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ABSTRACT 

We explore the thermal evolution of a neutron star undergoing episodes of intense accretion, separated by long 
periods of quiescence. By using an exact cooling code we follow in detail the flow of heat in the star due to the 
time-dependent accretion-induced heating from pycno-nuclear reactions in the stellar crust, to the surface photon 
emission, and the neutrino cooling. These models allow us to study the neutron stars of the Soft X-Ray Transients. 

In agreement with Brown, Bildsten and Rutledge (1998) we conclude that the soft component of the quiescent 
luminosity of Aql X-l, 4U 1608-522, and of the recently discovered SAX J1808.4, can be understood as thermal 
emission from a cooling neutron star with negligible neutrino emission. However, we show that in the case of 
Cen X-4, despite its long recurrence time, strong neutrino emission from the neutron star inner core is necessary 
to understand the observed low ratio of quiescent to outburst luminosity. This result implies that the neutron star 
in Cen X-4 is heavier than the one in the other systems and the pairing critical temperature T c in its center must be 
low enough (well below 10 9 K) to avoid a strong suppression of the neutrino emission. 

Subject headings: stars: neutron — X-rays: stars — dense matter — stars: individual (Aquila X- 1 , Cen X-4, 
SAX J1808.4-3658, 4U 1608-522, Rapid Burster, EXO 0748-676) 



1. INTRODUCTION 

Neutron stars in Soft X-Ray Transients (SXRTs) undergo re- 
current surges of activity separated by long phases of relative 
"quiescence". The detection of type I bursts during the X-ray 
brightening observed in outburst (~ 10 37 ~ 38 ergs _1 ) unambigu- 
ously indicate that episodes of intense accretion occur onto the 
stellar surface. The origin of the faint X-ray emission, ob- 
served in quiescence, at the level of ~ 10 32_33 ergs _1 , is instead 
uncertain and may result either from surface thermal emission 
(Brown et al. 1998) and/or from shock emission due to an hid- 
den fastly spinning pulsar (Stella et al. 1994; Campana et al. 
1998a,b). 

Brown, Bildsten and Rutledge (1998, BBR98) showed, using 
simple but concrete arguments, that the "rock-bottom" emis- 
sion, in quiescence, could result from the cooling of the neu- 
tron star made hot during the events of intense accretion: the 
inner crust compressed by the loaded material becomes the site 
of pycno-nuclear reactions that may deposit enough heat into 
the core to establish a thermal luminosity 
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detectable when accretion halts ( (M) is the time averaged ac- 
cretion rate, Q nuc ~ 1.5 MeV, the nuclear energy deposited per 
baryon, and m u the atomic mass unit). Following this consid- 
eration, Rutledge et al. (1999, 2000) showed that spectral fits 
with accurate hydrogen atmosphere models lead to emitting ar- 
eas consistent with a neutron star surface, pointing in favor of a 
thermal origin of the quiescent emission (or of part of it). 

The interpretation of the quiescent emission seen in SXRTs 
as due to "cooling" gives the opportunity of probing the physics 
of the neutron star interior in an unprecedented way. We study 
in detail the thermal evolution of transiently accreting sources 
with a full cooling code, computing the "exact" quiescent lumi- 



nosity and relating it to the accretion history. We then explore 
the consequences of the cooling hypothesis to infer properties 
of the underlying neutron star. 



2. THE COOLING MODEL 

We use an "exact" cooling code which solves the equations 
of heat transport and energy conservation in a wholly general 
relativistic scheme (Page 1989). The cooling sources are neu- 
trino emission in both the crust (electron-ion bremsstrahlung, 
plasma and Cooper pair formation) and the core (modified Urea 
and neutron/proton bremsstrahlung, Cooper pair formation, and 
direct Urea in some cases), and surface photon emission. The 
heating source is the accretion-induced production of nuclear 
energy. The energy released in the upper atmosphere by the in- 
falling matter and in the envelope by thermonuclear burning is 
here neglected since it is rapidly radiated away (Fujimoto et al 
1987; Brown 2000). The heat sources located in the crust are 
instead important (BBR98) and describe the readjustment, to 
chemical equilibrium, of the matter processed during the heavy 
accretion events. Energy release by electron captures and neu- 
tron emissions is included following Haensel & Zdunik (1990a) 
and proceeds at a rate proportional to the instantaneous accre- 
tion rate M(t ). Pycno-nuclear fusions, releasing the bulk of the 
energy, continue when the star is returning to its cooling state, 
since their characteristic time-scale (~ months) is comparable 
to or longer than the duration of a typical outburst. This time 
lag is accounted for with a suitable heat release function (see 
Possenti et al. 2001). The inclusion of all the non-equilibrium 
energy sources (down to p = 1.137 x 10 13 g cm" 3 ) assumes that 
the entire stellar crust has been replaced; this would take ~ 50 
million years for a time averaged accretion rate of 10 _1 l M® yr" 1 
(Haensel & Zdunik 1990a). The temperature gradients in the in- 
terior determine how much of this heat is stored and how much 
directly flows away to the surface. We calculate the thermal 
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conductivity using the results of Baiko and Yakovlev (1995) in 
the crystallized crustal layers, Itoh et al. (1983), Mitake et al 
(1984), and Yakovlev (1987) in the liquid ones, and Flowers & 
Itoh (1980) in the core. The calculations by Potekhin, Chabrier 
and Yakovlev (1997) of the structure of envelopes made of ac- 
creted material are used to relate the interior temperature, at a 
density of 10 10 g cm -3 , to the surface thermal luminosity L q , 
(i.e. the quiescent luminosity visible in between outbursts) and 
to the corresponding effective temperature T e s, 
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Fig. 1. — Critical temperatures T c for both proton 'So and neutron 'So and 
3 P2 pairing assumed in this work. The central densities of the stars of vari- 
ous mass we study are indicated as well as the threshold for the direct Urea 
("DURCA") process. 

The chemical composition and the equation of state (EOS) of 
the crust are the ones of an accreted crust (Haensel & Zdunik 
1990b) and the state of this matter as well as its specific heat 
are followed during the evolution of the star from the liquid 
to the solid phase (Slattery et al 1980) and the Debye regime 
(Shapiro & Teukolsky 1983). For the core we follow Prakash, 
Ainsworth & Lattimer (1988) using a compression modulus 
Kq = 240 MeV and a symmetry energy dependence oc p° 7 . With 
this particular EOS, the direct Urea process is allowed at densi- 
ties p > 1.28 x 10 15 g cm" 3 , which are attained in neutron stars 
heavier than a critical mass M cl - = 1 .65 M . 

Finally, we include the strong suppressing effects of super- 
fluidity (neutron) and superconductivity (protons) on both the 
specific heat and neutrino emission: given the enormous uncer- 
tainty (Baldo et al 1998) on the value of the critical temperature 
T c for 3 P2 neutron pairing we plot explicitly in Fig. [j] the values 
we adopt here. We will refer to our 1.4 and 1.6 M stars as 
slow cooling neutron stars and the 1.7 and 1.8 M Q ones as fast 
cooling stars: the former ones have only the modified Urea pro- 
cess allowed in their core and suppression by neutron pairing is 
strong while the latter ones have the direct Urea process operat- 
ing with suppression at lower temperatures so that fast neutrino 
cooling does affect their thermal evolution (Page 1998). 

3. TOWARD THE EQUILIBRIUM STATE 

As a preliminary step we verify whether transient accretion 
on an initially cold neutron star is able to bring it to a hot sta- 



tionary state and how fast this happens. 
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Fig. 2. — Redshifted core temperature T corc (10 8 K) vs time t soon after the 
onset of transient accretion, set at t = 0. The right scale reports the effective 
surface temperature as measured at infinity. The solid lines refer to a 1.4 Mq 
superfluid star with transient accretion, AM = 6 X lCr"M0 , ? rec = 150 days, 
and foui = 30 days. The upper curve shows the evolution from a hotter state that 
may result from an early phase of steady accretion, whereas the lower curve 
mimics a resurrecting accretion period following a phase of pure cooling. The 
thickness of the lines results from the rapid variations of the temperature due to 
accretion, as can be seen in the insert where T com is plotted against time over 
a few cycles. The dot-dashed line shows the evolution of the same star with 
continuous accretion at (M) = 1 .46 X 10~ 10 Mq yr -1 . 

For this, we model the accretion rate M(f) (as measured at 
infinity) with a fast exponential rise, on a time scale f r i se , reach- 
ing a maximum M max followed by a power law decay of index 
a = 3. Accretion is never turned off but becomes rapidly negli- 
gible in the decay phase. A new exponential rise recurs every 
free- frfse is set by the duration of the outburst phase f out (= 3f r i S e) 
which, for SXRTs is ~30 days long. The total accreted mass 
during one cycle, AM = 1 .08 f r i se M max , is the key parameter and 
is used to determine M max . We define the outburst luminosity 
(at infinity) L out such that the fluence is L out f out = r\ AMc 2 with 
an efficiency ?7 = 1-e^ (e^ = (1 -IGM/Rc 1 ) 1 ! 2 ~ 0.15-0.30 
is the redshift factor). The time averaged accretion rate is 
(M) = AM/f rec . 

We show in Fig. ^ an example of the heating (cooling) of an 
initially cold (hot) neutron star after the onset of intermittent ac- 
cretion. It is compared with a model having constant accretion 
at a rate equal to (M). As expected, the star reaches thermal 
equilibrium and does it on a very short time scale r equ ~ 10 4 yr, 
much shorter than any binary evolution time, and, in particu- 
lar, much shorter than the time necessary to replenish the crust 
with fresh non-catalyzed matter. The equilibration temperature 
is reached when the net injected heat is exactly balanced by the 
energy loss from the surface and/or from neutrino emission. 

4. THE QUIESCENT LUMINOSITY 

Our stars reach a hot equilibrium state and glow, in quies- 
cence, at a luminosity L q which varies with f lec and AM, as 
illustrated in Fig. |3|. The solid lines refer to models with AM = 
6 x 10 _I1 M Q , compatible with that inferred for Aql X-l, 4U 
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1608-52, and Cen X-4. For a 1.4 Mq star we have also drawn 
the curve corresponding to a lower value of AM = 1O -11 M0 , 
that may describe the fainter transient SAX J 1808.4. L q de- 
pends crucially on whether fast neutrino emission in the inner 
core is allowed or inhibited, as illustrated in Fig. [|. 

For slow cooling stars with superfluidity (upper solid lines 
'1.4sf ' and '1.6sf'), neutrino emission is totally suppressed and 
the luminosity is determined by balance of nuclear heating with 
photon cooling. When superfluidity is not included, these stars 
are slightly less luminous as there is a small leakage of neutri- 
nos due to the inefficient modified Urea process; the difference 
is very small and we omit the results in the figure. 

When fast neutrino emission is allowed ('1.7sf and '1.8sf 
lines), the luminosity is much lower and the lowest L q are ob- 
tained in the case neutrino emission is not affected by neutron 
pairing ('l.Snsf case). In these fast cooling models the depen- 
dence of L q on f rec is much weaker than in the slow cooling 
cases since most of the heat deposited during an accretion cy- 
cle is rapidly lost into neutrinos. If T c were very high at the 
center of even our most massive stars, then neutrino emission 
would be completely suppressed and the evolution of the su- 
perfluid 1.8 Mq star would be almost identical to the super- 
fluid 1 .4 M Q star. This is similar to what happens in cooling 
isolated neutron stars: the star's temperature is almost entirely 
controlled by the value of T c (Page & Applegate 1992). We 
explicitly choose a relatively low T c at high density to allow 
heavier stars (1.7 and 1.8 Mq in our specific model) to have a 
qualitatively different behavior relative to the less massive ones. 
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Fig. 3. — The (equilibrium) quiescent luminosity versus the recurrence time 
resulting from our numerical calculations. The dashed area covers the region 
of the observed luminosity in SXRTs. The solid lines refer to a star accreting 
AA/ = 6 X 10 _1 'Mq in each cycle. They are labelled by the mass of the star (in 
solar masses) and the existence (sf) or nonexistence (nsf) of a superfluid phase. 
The dot-dashed line shows our results for a 1.4 Mq superfluid star loading 
AM = 10 _1 1 Mq in each cycle. 



5. COMPARISON WITH DATA 

The most convenient quantity for comparison with data is 
the ratio 1Z of fluence in quiescence L q f rec to fluence in outburst 



^out ^out 

~ rjAMc 2 (BBR98). This gives us a direct measure- 
ment of the fraction / of the heat released in the crust during 
accretion which, stored in the stellar interior, slowly leaks out 
to the surface 
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Our exact cooling code computes L q and thus predicts /. This 
is illustrated in Fig. Q where we give model values for L q /L out 
against t lec /t out and where we draw the theoretical line /= 1 
of Eq. |2] for 1.4Mq and I.8M0 stars. Our numerical results 
when neutrino emission is negligible ("1.4sf" solid line) are 
very close to the / = 1 line showing a heat storage efficiency 
almost independent of t rec /t out . The ratio L q /L out is dramati- 
cally reduced when direct Urea emission switches on, as in the 
core of our heavier stars ("1.7sf" and "1.8sf" solid lines). In 
these cases, the heat storage efficiency decreases with decreas- 
ing ?recAout (along a line of given mass): the neutron star in a 
transient with shorter f rec has a warmer core and an enhanced 
neutrino emissivity (oc T^ OK ) leading to a lower value of /. 
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Fig. 4. — Quiescent to outburst luminosity ratio plotted versus the ratio of 
the recurrence time over the outburst time. Bold solid lines are the results of 
our numerical calculations as in FigM The expected results for / = 1, see Eq.bl 
are shown for our 1 .4 Mq star (dashed-dot line) and 1 . 8 Mq star (dashed line). 
The filled squares represent the observed values. 

Rutledge et al. (2000) analyzed the observation of five 
SXRTs, Aql X-l, Cen X-4, 4U 1608-52, EXO 0748-767, and 
the Rapid Burster, made during the periods of lowest emis- 
sion. Fitting the spectra with the H atmospheric models of 
Zavlin et al (1996) they calculated the thermal luminosity L q 
and the outburst time averaged luminosity L out . All five sys- 
tems display similar peak luminosities and outburst times. For 
the fainter source SAX J1808.4, we take L q = 2.5 x 10 32 ergs _I , 

^out — 

erg s , f ou t = 20 days and ?, ec = 700 days (Stella et al. 
2000). Overlying all these observations in the diagram of Fig. [|, 
it appears remarkable that Aql X-l and 4U 1608-52 stands just 
on the theoretical line constructed for 1 .4 Mq (close also to 1 .6 
Mq as shown in Fig. The Rapid Burster would fit the theory 
if its mass lies in an intermediate range with fast neutrino emis- 
sion strongly controlled by pairing. Cen X-4 is adjacent to the 
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1.8 Mq model star results, having an /i.g of 0.18. Such a low 
efficiency means that about 80% of the deposited heat is not ob- 
served to flow back to the stellar surface: only a very efficient 
UnRecorded Cooling Agent such as the direct Urea process is 
able to provide such a strong energy drain. Note however that 
the ratio t rec /t out , for Cen X-4, could be underestimated as only 
two outbursts have been recorded, in 1969 and 1979. Thus, we 
can not exclude that Cen X-4 is close to the curve /1.4 = 1; this 
however requires a very large t TSX /t ont of 1,000. Its low value of 
L q /L out could alternatively indicate the presence of a not fully 
replaced crust. We here take the view that Cen X-4 is heavy. 

Having calculated the upper theoretical bound on the effi- 
ciency of the rediffusion of the pycno-nuclear energy, we can 
interpret the peculiar behavior of EXO 0748-767 as due to 
an extra luminosity resulting from the gravitational energy re- 
leased by a faint accretion and/or by the interaction of the in- 
falling matter with the magnetosphere (BBR98). 

Our results confirm the overall picture suggested by BBR98 
in their parametric approach: the quiescent luminosity seen in 
SXRTs comes from the rediffusion toward the surface of the 
heat deposited in the core by the pycno-nuclear reactions trig- 
gered in the crust due to transient accretion. 

6. WEIGHTING NEUTRON STARS AND PROBING SUPERFLUIDITY 

We can now use the established link of / with mass, to infer 
properties of the star in a transient through the comparison with 
the data. The mass signature in the L q /L out vs t Itx /t out plot is 
mainly determined by the state of the matter in the inner core 
of the neutron star. The ratio L q /L out has a clear drop when the 
source mass is above the critical mass for fast neutrino emis- 
sion M CI determined by the value of the pairing critical temper- 
ature T c in the innermost core. The value of M cr =1.65 Mq 
used here is model dependent, but combining future indepen- 



dent measurements of the mass with the observed ratio L q /L out 
will open the possibility of actually measuring the value of M cr 
and/or also infer an upper limit on the critical pairing tempera- 
ture r c . Notice, however, that many other fast neutrino emission 
channels, beside the simplest nucleon direct Urea considered 
here, are possible and the constraints we could obtain about the 
value of T c may then well refer to hyperon or quark pairing 
(Page et al. 2000). 

Awaiting future independent mass measurements, we can 
discriminate between stars with masses close to the canonical 
value 1.35-1.40 M Q (Thorsett & Chakrabarty 1999) and more 
massive objects. In this context, the ratio L q /L out of Cen X-4 
can be explained if its neutron star is heavier, and even probably 
much heavier, than in the other observed SXRTs. In low mass 
binaries, the neutron star could have undergone a substantial 
mass load (0.1 to 0.5 M Q ) while being spun up to very short 
periods. Observational evidence (Casares et al. 1998) and the- 
oretical arguments (Stella & Vietri 1999; Possenti et al. 1999) 
can support this view. 

Charting the temperature of the old hot neutron star in a soft 
X-ray transient is hence a valuable tool to investigate dense 
matter in neutron star cores which is complementary to the 
study of isolated young cooling neutron star. It moreover opens 
the possibility to study the interior of stars which can be more 
massive than the isolated ones. 
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